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Complex systems may be subject to various uncertainties. A great 
effort has been concentrated on predicting the dynamics under uncertainty 
in initial conditions. In the present work, we consider the well-known 
Burgers equation with random boundary forcing or with random body 
forcing. Our goal is to attempt to understand the stochastic Burgers 
dynamics by predicting or estimating the solution processes in various 
diagnostic metrics, such as mean length scale, correlation function and 
mean energy. First, for the linearized model, we observe that the important 
statistical quantities like mean energy or correlation functions are the same 
for the two types of random forcing, even though the solutions behave very 
differently. Second, for the full nonlinear model, we estimate the mean 
energy for various types of random body forcing, highlighting the different 
impact on the overall dynamics of space-time white noises, trace class 
white-in-time and colored-in-space noises, point noises, additive noises or 
multiplicative noises. 
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1 Introduction 



The Burgers equation has been used as a simplified prototype model for hy- 
drodynamics and infinite dimensional systems. It is often regarded as a one- 
dimensional Navier-Stokes equation. Our motivation for considering this equa- 
tion comes from the modeling of the hydrodynamics and thermodynamics of 
the coupled atmosphere-ocean system. At the air-sea interface, the atmosphere 
and ocean interact through heat flux and freshwater flux with a fair amount 
of uncertainty |361 1191 IHUj . These translate into random Neumann boundary 
conditions for temperature or salinity. The Dirichlet boundary condition is also 
appropriate under other physical situations. The fluctuating wind stress forcing 
corresponds to a random body forcing for the fluid velocity field. The coupled 
atmosphere-ocean system is quite complicated and numerical simulation is the 
usual approach at this time. In this paper, we consider a simplified model 
for this system, i.e., we consider the Burgers equation with random Neumann 
boundary conditions and random body forcing. Although the stochastic Burg- 
ers equation is widely studied, most work we know are for Dirichlet boundary 
conditions or periodic boundary conditions |181 1231 ITT1 112j . The reference |37| 
studied the control of deterministic Burgers equation with Neumann boundary 
conditions. 

We consider the stochastic Burgers equation with boundary forcing on the 
interval [0, L] 

dtu + u ■ d x u = vd\u (1) 

a x u(-,0)=a77 a x u(-,L) = 0. (2) 

Here a > denotes the noise strength and r] is white noise, i.e., rj is a generalized 
Gaussian process with Er/(i) = and Krj(t)rj(s) = 5{t — s). The restriction to 
noise on the left boundary is only for simplicity. Analogous results will be true, 
if forces act on both sides of the domain. 

We will see that boundary forcing coincides with point forcing at the bound- 
ary. Thus we also look at point forcing. As a simple example for point forcing, 
we consider 

dtv + v ■ d x v = vd^v + a.5orj (3) 

U(;-L)=U(;L)=0, (4) 

where 5o is the Delta-distribution. 

We will compare solutions of (0) and © with solutions of the stochastic 
Burgers equation with body forcing. 

dtv + v ■ d x v = vd 2 £ v + <r£ (5) 

either subject to Dirichlet or Neumann boundary conditions. Here the noise 
strength is denoted by o > and £ is space-time white noise. I.e., £ is a 
generalized Gaussian process with E£(i, x) = and E£(£, x)£(s, y) = S(t — 
s)5(x — y). We will also consider trace class body noise, i.e., noise that is white 
in time but colored in space. 
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For the linearized equations, we will compare statistical quantities of both 
solutions, which are frequently used. One of them is the mean energy 

1 l' L 

- \ E[u{t,x) -u(t)] 2 dx, (6) 
L Jo 

where 

_ 1 f L 

u(t) = — I u(t, x)dx. 

L Jo 

Another important quantity, which gives information about the characteristic 
size of pattern, is the mean correlation function 



1 



L 



C(t,r) := - E[u(t,x)-u(t)]-[u(t,x + r)-u(t)]dx, (7) 
L Jo 

which is usually averaged over all points r with a given distance from 0. We 
obtain the averaged mean correlation function 

C(t,r) = ±[C(t,r) + C(t,-r)] (8) 

where we employ the canonical odd and 2L-periodic extension of u in order to 
define C(t,r) for any r£l. 

For the linearized equation the main result states that mean energy and 
averaged mean correlation function are the same for solutions of Q and Q. 
Nevertheless the solutions behave completely different. Furthermore, we give 
some qualitative properties like, for instance, the typical pattern size. This 
should carry over to a transient regime (i.e., small times) for the corresponding 
nonlinear equations. 

For the full nonlinear Burgers model, we estimate the mean energy for vari- 
ous types of random body forcing, highlighting the different impact on the over- 
all dynamics of space-time white noises, trace class white-in-time and colored- 
in-space noises, point noises, additive noises or multiplicative noises. 

In the following, we discuss linear dynamics in §2 and nonlinear dynamics 
in S3. 



2 Linear Theory 

Define 

A = ud 2 x 

with 

D(A) = {w£ H 2 ([0,L}) : 8 x w(0) = 0, d x w(L) = 0} 

It is well-known (cf. e.g. ^H]) that A has an orthonormal base of eigenfunctions 
{efcjfceNo i n L 2 ([Q., L\) with corresponding eigenvalues {Afc}fc e rsj . In our situation 
eo(x) = 1/V~L, efe(x) = y/2/L ■ cos(irkx/L) for k & N, and Afc = — is(k7r/L) 2 . 
Moreover A generates an analytic semigroup {e tA }t>o- (cf. e.g. |3*T]). 
In fact, e tA vo is the solution of the following evolution problem 

d t v = Av, d x v(-,0) = d x v(; L) = 0, u(0, x) = v (x). (9) 
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The solution is 

oo 

e tA v (x) := v(t, x) = V] < v Q , e k > e Xkt e k , t > 0, < x < L, (10) 

k=0 

where < •, • > is the usual scalar product in L 2 (0, L). 

We now consider the following linearized problems. First 

d t u = Au, d x u(-, 0) = ad t /3, d x u(-,L)=0. (11) 

Here the white noise rj is given by the generalized derivative of a standard 
Brownian motion (cf. e.g. [Q), and a is the noise intensity. 
Secondly, 

d t v = Av + ad t W, d x v(-,0) = 0, d x v(-,L) = 0, (12) 

where the space-time white noise is given by the generalized derivative of an 
Id- Wiener process. Namely, W(t) = X^fceN h{t) e ki where {(3 k }keN is a family 
of independent standard Brownian motions, and a is the noise intensity. 

It is known (cf. e.g. ^1) that (|12j) has a unique weak solution given by the 
stochastic convolution (taking initial condition to be zero) 

W A (t) =a- I e^ A dW(r) = a ■ V / df3 k {r)e k . (13) 

We define the Neumann map T> by 

(1 - A)V~f = 0, d x P 7 (0) = 7 , d x V 7 (L) = 

for any 7 G R. It is known that V : R i-> H 2 ([0,L]) is a continuous linear 
operator. In fact, we have explicit expression for this linear operator 

pX _i_ p2Z/ p x 

ni) = l_J L 7. (14) 

From ^B] or we immediately obtain, that (|11|) has a unique weak solution 
(taking initial condition to be zero) 

Z(t) = (l-A) I ^-^ A Vad{3(T). (15) 
J 

In the next section we derive explicit formulas for Z in term of Fourier series. 
2.1 Mean Energy 

To obtain the Fourier series expansion for Z, consider for e G D(A) and 7 G R 

< 2?( 7 ), (1 - A)e >£2( [0 ,£]) = <2?( 7 ),e>- / V( 7 ) ■ e xx dx 

Jo 

f L 

= <V{i),e>-l T>(j) xx ■ edx + V{i) x ■ e|^~ 

JO 

= - 7 e(0), (16) 
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by the definition of T>. Hence, 



<Z(t),e k > = < [ e^ A Vad/3(T),(l-A)e k > 
Jo 



< Vad/3(r),(l - A)e k > 



ae k (0) 



Jt-r)X k 



dp(r). 



We now obtain 



Z(t) = a- J> fc (0) f 



d/3(r)e k . 



Finally, 



and 



Z(t) = aei(O) • [ e^ T)Xk df3{T)e k + aeg(0)/3(t) 
W A (t) = a • V f e^^dp k {r)e k + a/3 (i) . 



If we now choose a = aei(O), we readily obtain that 



(17) 



(18) 



(19) 



(20) 



E\\Z(t) - Z(t)\\ 2 = <J 2 -J2 te 2rXk dT = E\\W A (t) - W A (t)\\ 2 , 

km 

where || • || is the norm in L 2 ([0,T]). Hence, the mean energy in both cases is 
given by a 2 ^ 1 £ fceN /„* e 2rX "dT. 

For the mean energy we can prove the following theorem, which is similar 
to the results of [I] and jSJ. 

Theorem 1 Fix a 2 = a 2 /L, then the mean energy Cz(t,0) = Cyy A (t,0) be- 
haves like Ci(a 2 /L)^t/u for t <C L 2 jv, and like C20? jv for t » L 2 /u. 

The main difference to body forcing is the scaling in the length-scale L. The 
long-time scaling is independent of L, while the transient scaling is. 

C z (t,0) 




Figure 1: The scaling of the mean energy for boundary forcing. 
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2.2 Correlation Function 

To obtain results for the correlation function, we think of Z and Wa to be 
periodic on [— L, L], and symmetric w.r.t. 0. I.e., we choose the standard 2L- 
periodic extension respecting the Neumann boundary conditions on [0, L\. To 
be more precise, we extend Z and Wa in a Fourier series in the basis e^, which 
we then consider to be defined on whole K. 
We consider firstly for k, I ^ 



L 

ek(x)ei(x + r)dx 



Le k (0)e t (0) ^fl^ sm^lr/L) : k + l 
Lefc(0) 2 cos(7rkr/L) : k = l 



Now relying on the independence of the Brownian motions, it is straightforward 
to verify 

Cw A {t,r) = ^■E<W A (t,x)-W A ~{t),W A (t,x + r)-W A -(t)> 



a 2 ej{0) 



L 



V f e 2rXk dr cos{itkr / L) , (21) 



as a 2 ef (0) = a 2 . Furthermore, 

C z {t,r) = — E < Z(t, x) — ~Z(t), Z(t, x + r) — Z(t) > (22) 
= C WA (t,r) + ^ £ /e^ ^V-^ 2 ) ^ S[n{7rlr/L) ' 

k^l 

Obviously, Cz and Cw A do not coincide, but let us now look at the averaged 
correlation function 

C(t,r) = ±[C(t,r) + C(t,-r)}. 

Then it is obvious that 

C WA (t,r) = C WA (t,r) = C z (t,r) ± C z {t,r) . (23) 

Now 

Theorem 2 For a 2 e 2 (0) = a 2 the mean energy and the averaged mean corre- 
lation functions Cw A and Cz for Z and Wa coincide for any t > 0. 

This is somewhat surprising, as realizations of Z and Wa behave completely 
different, when the condition a 2 ef(0) = a 2 is satisfied. See e.g. Figure HI and 
Figure 01 

It is even more surprising, as the scaling behavior of quantities like mean 
energy and mean correlation functions are an important tool in applied science, 
which for example is used to determine the size of characteristic length scales 
and the universality class the model belongs to. Here both linear models lie in 
the same class, although their behavior differs completely. 
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Figure 2: Random boundary condition: One realization of the solution of the 
equation for L = 1, v = l, a = l and initial condition u(x,0) = 0. 



The scaling behavior with respect to L and t of the mean energy can be 
described using the results of [3] , where the mean surface width for very general 
models was discussed. Therefore we focus on the scaling properties of the mean 
correlation function. Here we also want to investigate the dependence on a and 
v. 

First we consider the scaling properties of the correlation function Cz(t,r) 
or Cw A , as given in (|2*2*)) . We are especially interested in the smallest zero of the 
function, which gives information about characteristic length scales or pattern 
sizes. For this, we use the normalized correlation function. 



Pz(t,r) 



Cz(t,r) 
C z (t,0) 



(24) 



Note that Cz{t, 0) is the mean energy and the maximum of Cz(t, r). 

We begin with some technical results. For any continuously differentiable 
and integrable function / : M + -tKwe obtain using the mean value theorem 



/ f(x)dx-J2f(k) 
Jo k=i 

For f(k) := e ~ 2Tuk2n2 1 1,2 cos(/c7rr/L) it is easy to verify that 



<Y, ^p \f'( v )\. 



(25) 



\f(k)\ < e -^ 2 - 2 ^ 2 . 1 



r + 



where we used that x s e x2 ° < (2ae) 1 I 2 for any x, a > 0. Hence 

4 



V sup |/'(t7)| 

fc=1 r?e(fc-i,*0 



oo 

k=l ^ 



'2e 



TV\ 
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Figure 3: Random body forcing: One realization of the solution of the equation 
(jl2() for L = 1, u = 1, (7 = 1 and initial condition v(x, 0) = 0. 



ir r 4 . ., f c 

£ i |r + 7S^ l(1 + X 

= -[r + -=Vru\{l + — p= 



-TUk 2 TT 2 /L 2 



dk) 



and 



Moreover, 



e~ fe dfc) (26) 

/ V sup |/'(»/)|dr<C T J-[r + v^][L + V^]. (27) 
Jo fr[ V &(k-l,k) l> v V 



t roo 



JO 



f{k)dkdr 



1 

7T 

1 /t 



1 - e 



-2tuk A 



cos(kr)dk 



o 



2isk 2 
r 

G{—=) 



with G(x) := 1 | fc2 — cos(kx)dk. 

1.6 
1.2 
0.8 
0.4 

1 2 3 4 5 6 

Figure 4: A sketch of G 
Using (|23j) we immediately obtain 




(28) 



Note that the approximation with G is not L-periodic in r, while Cz(t,r) is. 
The solution is that the error term is 0(1) for r near L. 
For the normalized correlation function we deduce 

, = C z (t,r) = G(^) + O (U- + Vt^][L + yg) 
PZ( ' rj ' C z (*,0) G(O)+0(^[v^][L + v^]) 

From the properties of G we infer the following: 

Theorem 3 Given 5 £ (0, 1) and sufficiently small €2 > 0, t/iere exists some 
ei > and three constants < Ci < C2 < C3 depending only on 5 and €2 such 
that for t < e\L 2 jv the following holds: 

pz(t,r)>8 for re[0,CiVH 

and 

\pz(t,r)\ < e 2 for r G [C 2 V^, C 3 V^]. 

Note that we did not show that the correlation function has a zero, but it is 
arbitrary small in a point r e rs Therefor the theorem says that the typical 
length-scale is y/tu, at least for times t <C L 2 /i/. 
For t — > 00 we immediately obtain that 



and 




Figure 5: A sketch of F 

We can look for the explicit representation of F, which is a 2-periodic func- 
tion, and compute explicitly the zero, but all we need from F is, that for a given 
small enough 5 > there is a x$ > such that F > 5 on [0,2:5]. Moreover, 
there is some xq such that F(xq) = 0. 
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Consider the normalized correlation function 



C z (t,r) F(r/L)+Q(e- 2t ™ 2 / L2 ) 
PZ[ ,r) C z (t,0) F(0) + O(e-^ 2 /i 2 ) ■ 

Assume that tv S> L 2 (i.e., there is some small e > such that etu > L 2 ). Now, 

pz(t,x L) = O(e- 2t ^ 2 / L2 ) 

and 

pz(t,xL)>^- )+ O(e~ 2t ^ 2 / L2 )>0 

for any x < x$. 

So for t S> L 2 /V the first zero of pz should be of order L. A more precise 
formulation is: 



Theorem 4 Given 5 G (0, 0.8) and 5 3> e2 > 0, i/iere existe some ei > 0, a 
constant C > 0, and a pomi x > depending only on 5 and €2 such that for 
t > L 2 j{ye\) we obtain the following: 

pz(t,r)>5 for r G [0, CL] 

and 

\p z (t,x L)\ < e 2 . 

Thus the theorem tells us that for t 3> L 2 /v, the typical length scale is of order 
L, which is the size of system. This result is true for both boundary and body 
forcing. 



3 Nonlinear theory 

For the nonlinear results we leave the setting of boundary forcing. Mainly, due 
to the lack of a-priori estimates. Usually, for Neumann boundary conditions 
one relies on the maximum principle to bound solutions, but the solution for 
boundary forcing is quite rough. Therefore, we hardly get useful results. Only, 
transient bounds for small times are possible to establish. For the next sections, 
we focus first on body forcing and later on point forcing. We will see later that 
boundary forcing is actually just a point forcing in a point at the boundary. 

The main results of this sections are uniform bounds on the energy and thus 
on the correlation function C, as |C(i,r)| < C(i,0), and C(t,0) is the energy. 
Furthermore, we show that for t — > the linear regime dominates. In also 
H" older-continuity for the mean energy was shown for a quasigeostrophic model. 
We conclude this section by a qualitative discussion on upper bounds for the 
energy using additive and multiplicative trace-class noise. 
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3.1 Body forcing - Mean energy bounds 

Here we provide bounds on the mean energy for the body forcing case. We 
consider additive space-time white noise case first, and show that the mean 
energy and thus the correlation function is uniformly bounded in time. This 
result is known (cf. [29 ) for Burgers equation using the celebrated Cole-Hopf 
transformation, but we provide here a simple proof for completeness. Further- 
more, our proof is based on energy estimates and it is easily adapted to other 
types of equations and additional terms in the equation. In contrast to that 
Cole-Hopf transformation is strictly limited to the standard Burgers equation. 

For a long time for space-time white noise only uniform bounds for log- 
arithmic moments were known. See |16| Lemma 14.4.1] or P3]. In [201 the 
transformation to a stochastic heat equation via the celebrated Cole-Hopf trans- 
formation was used to study finiteness of moments. Here we rely on a much 
simpler tool, which can also be applied to other equations. See for instance [Hj 
for a quasigeostrophic model, where our analysis would apply, too. 

Consider 

d t u + u ■ d x u = vdlu + ad t W , (29) 

u(-,-L) =u(-,L) =0, u(x,Q) =u (x). (30) 

Here W is a Q- Wiener process with a continuous operator Q G £{L 2 ). Thus 
W might be cylindrical, and we include the case of space-time white noise. 
Using the semigroup e l A the solution for this system is (see |16| or P.8 ): 

u(t) = e tA u - f e {t ~ T)A (\$x(r) + ld x u(r,x) 2 )dT + $ x (t) , (31) 

where for some A > fixed later 

<S> x (t) = aa f e^ T ^ A -^dW{T) 
Jo 

solves 

d t $ = vdl® - A$ + ad t W 

subject to 

H-,-L) = $(-,L) = 0, $(x,0)=0. 

Our main result is now: 

Theorem 5 Consider initial conditions uo with E\\uo\\ 2 < oo, which are inde- 
pendent of the Wiener process W (e.g. deterministic). Then the mean energy 
of the solution of hHl\) is uniformly bounded in time. I.e., 

supE||w(i) - ?l(t)|| 2 < oo . 
t>o 

Remark 1 Actually, we prove that sup i>0 E||ti(t)|| 2 < oo. The main problem 
in the proof is that after applying Gronwall-type estimates we end up with terms 
Eexp{ J * || < I > a(s)|||oo}. This might blow up in finite time, as second order expo- 
nential moments of the Gaussian may fail to exist, if t is too large. This is why 
we introduced artificially additional dissipation in the equation for &\, in order 
to get exponential moments small. 
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For the proof of Theorem [5] define 

v(t) = u(t) - * A (t) for t > 0, A > . (32) 
We see that v is a weak solution of 

d t v + ^d x (v + <Z> x ) 2 = vd 2 x v + \$> x (33) 

v(-,-L) = v(-,L)=0, v(x,0)=u {x). (34) 

The following calculation is now only formal, but it can easily made rigorous 
using for instance spectral Galerkin approximations. Taking the scalar product 
in (|3"3*|) yields 



-dt\H 2 
2 11 11 



/L rL 
(v + <S>\fv x dx+ I <& x v dx 
L J-L 

< -\\v x \\ 2 + H^aII^-iIKI! + I^||||^a|||4 + 2||« a .[|[|$Alk o °[|« I 



li 4 

< ~c 2 p \\v\\ 2 + 4|M| 2 ||<I> A || 2 + 2\ 2 \\<S> X \\ 4 L4 + 211^11^ , 

where we used Young inequality (ab < |a 2 + 2b 2 ), and Poincare-inequality 
IMI < CpH^ || . Now, from Gronwall-type inequalities 

IK*)|| 2 < e- c p t+8 /oll*A|[L^|| u ( )||2 ( 35 ) 

JO 

Now we use the following lemma, which is easily proved by Fernique's theorem, 
if we consider <3?a 

Gaussian in L 2 ([0, t ], L°°). 
Lemma 1 Fix K > and to > 0, then there is a Xq such that 

sup Eexp{16 f \\$x(s)\\ 2 Loo ds} < K 2 
te[o,t ] J o 

for all A > Ao- 

Furthermore, we use that all moments of ||$a||l°° an d H^aIIh- 1 are uni- 
formly bounded in time. This is easily proven, using for instance the celebrated 
factorization method. 

Now we first fix K > 0, and then to such that e~ c p t K < \. This yields for 
t G [0, to] and A sufficiently large 

IE|Kt)|| 2 < e- c ^E||n(0)|| 2 +4K^- c p(*- s )(E(A 2 ||$A|ll4 + ||$A|li-i) 2 ) 1/2 ds , 

using Holder, LemmaEI and the independence of u(0) from <I>a. We now find a 
constant C depending on to and K such that 

sup E||t>(t)|| 2 < i^E||u(0)|| 2 + C and E||w(t )|| 2 < ^E||n(0)|| 2 + C 

te[o,t ] 4 
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Using E||ii(t)|| 2 < 2E||v(t)|| 2 + 2E||$ A (i)|| 2 yields for a different constant C 

sup E\\u(t)\\ 2 < KE\\u(0)\\ 2 + C and E||u(t )|| 2 < ^E||u(0)|| 2 + C 
te[o,t ] 2 

Now we repeat the argument for k G N by defining v(t) = u(kto + t) — $\(t), 
where $ A (i) has the same distribution than <J? A (i) due to a time shift of the 
Brownian motion. Now v solves again (|33|) with initial condition u{kto). Note 
that by construction u(kto) is independent of <£ A . 

Repeating the arguments as before yields for k £ No 

sup E\\u(t + kt )\\ 2 < KE\\u(kt )\\ 2 + C 
*e[o,to] 

and 

E|| U ((A; + l)to)|| 2 < ^E\\u(kt )\\ 2 + C . 

Now the following lemma, which is a trivial statement on discrete dynamical 
systems, finishes the proof. 

Lemma 2 Suppose for q < 1 and some C > we have a n+ \ < qa n + C , then 
a n is bounded by 

C 

a n < h a . 

1-q 

3.2 Point forcing - Mean energy bounds 

Consider hyperviscous Burgers equation with point-forcing. We would like to 
proceed exactly the way, we did in the previous section, But we can not, as for 
point forcing, the solution of the linear equation might fail to be in L°°. This 
is why we add additional damping. Hyperviscous Burgers equation has been 
studied in several occasions. See for example [71 [2*71 152*] . 

Consider for some e > the operator A e = — v{— A) 1+e , where A is the 
Laplacian subject to Dirichlet boundary conditions. Then the hyperviscous 
Burgers equation is given by 

dtu + u ■ d x u = A e u + a5o[3 (36) 

u(-,-L) =u(-,L) = 0, u(x,0) = it (x). (37) 

Here, (3 is a standard Brownian motion and 5q the Delta-distribution. 

Using the semigroup e tAt the solution for this system is (see [Ej or |18j): 

u(t) = e tA *u - f e^- T ^ (A* A (r) + \d x u(r, x) 2 )dr + $ A (t) (38) 
Jo 1 

where for some A > fixed later 

* A (*) = a f e (t - r ^- X h (x)dP(r) 
Jo 
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solves 

d t <& = A e $ - A$ + a5 $ 

subject to 

$(.,-L) = *(.,L) = 0, <J>(x,0)=0. 



Using the standard orthonormal basis {ek}keN of eigenfunctions of A e given by 

2L 



efc(x) = -y/l/L-sin(— with corresponding eigenvalues Afc = — (7rfc/2L) 2+2e , 



we see 



<D A (t) = aV / e(*" T 



-r)(A fc -A) d/3(r)efc(0)efc _ (39) 



Note that the Fourier-coefficients of that series are not at all independent. 
Thus we cannot rely on the better regularity results available for the stochastic 
convolution of the previous chapter. Especially, for e = we cannot show that 
*x(t)eL°°([-L,L]). 

Note that the series expansion of boundary and point forcing is very similar. 
Thus we can regard boundary forcing at a point forcing at the boundary, when 
the equation is subject to Neumann boundary conditions. 

Our main result is now: 

Theorem 6 For all e > and all initial conditions uq independent of (5 with 
i?||uo|| 2 < oo the solution of Ofll) satisfies that the mean energy is uniformly 
bounded in time. I.e., 

supE||ti(i) -u(t)\\ 2 < oo . 

We will proceed exactly as in the previous section. Now v = u — <3?a is a weak 
solution of 

d t v + ^d x (v + <S> x ) 2 = vd 2 x v + \$ x , (40) 

again subject to Dirichlet boundary conditions and initial condition v (0) = no- 
Now consider first the nonlinear term for some small 5 > 0. Using Holder, 
Sobolev embedding of H 2 p into LP and the bound 

||u||H*r < C\\A2 /(1+e) u\\ 

yields 



L 

v^xVxdx < ||u||i2+a||$A|li(4+25)/«||w a 



1 1 <5 11 
>A 1 1 1/4+2,5)/* | 



< c\\Ar + ' 2+s v\m x \\ T u +2 s )/s \\Ar +e '. 



Now we can easily find an 5 > sufficiently small such that there is a p 
p(e) £ (2, oo) such that (using interpolation inequality) 

L I 
v^xvxdx < C\\v\\z?\\$\\\Lp\\Aiv\\ . 
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Now we can use the same proof as in the section before. We only need that 
&\(t) £ L°°(0, L). To be more precise, an easy calculation using the series 
expansion of (|39|) shows that for any e > 

supE||$ A (£)||!oo < C e supE||$ A (t)|| 2 i +6 -> as A -> oo . 

t>0 ' t>0 H~2~ 

It is now straightforward to prove an analog to Lemma ^ The remainder of 
the proof is analogous to the section before. 

Let us remark that we could even simplify that proof a little bit, by avoiding 
second order exponentials of $ A . In that case we could work with A = 

3.3 Body forcing - Transient Behavior 

Let us focus on Burgers equation with body forcing. The results for hyper- 
viscous Burgers with point-forcing are completely analogous. We will prove: 

Theorem 7 Let u be a solution of \29\l and consider for simplicity u(0) = 0. 
Denote by 

E u {t) =E\\u(t)-u{t)\\ 2 
the mean energy ofu(t), then there is some So such that 

E u {t) = E$ (t) + 0(t3+*>) fort^O. 

To be more precise, for some to > sufficiently small there is a constant C > 
such that \E u (t) - | < Ct^ +S ° for all t G [0,i ]- 

As we know from results like Theorem ^ that E$ (t) behaves like y/t for small 
t, we can conclude that the linear regime dominates for small t. 

We could explicitly calculate So, but omit this for simplicity of presentation. 

For the proof of Theorem [7| use 

- E 9o (t)\ <CE\\v(t)\\ 2 , 

where we used Cauchy-Schwarz inequality and uniform bounds on E||n(t)|| 2 and 
E||<3? (i)|| 2 . Using (j3SJ) with A = and u(0) =0 yields together with Lemma [T] 

E|Kt)|| 2 <C /"(E||*o(t)|&-0 1/2 *- 
Jo 

It is now easy to show that E||3>o(*)||jj-i behaves like t 2S ° for some Sq > 0, 
which can be explicitly calculated using the methods of Theorem ^ Theorem 
|7| is now proved. 

A simple corollary using Holders inequality is: 

Corollary 1 Under the assumptions of Theorem^ we know for the mean cor- 
relation function 

C u (t, r) = C$ (t, r) + O(t^ +So ) for t^O and all r . 

Notice that this result is only useful for small times and small r, as seen from 
the qualitative behavior of C$ , which is similar to the results shown in section 
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3.4 Trace class noise: Additive vs. multiplicative body noises 

Consider again a solution of the following Burgers equation: 

d t u + u ■ d x u = ud 2 u + aW (41) 

u(; 0) = 0, u(; L) = 0, u(x, 0) = uq(x), (42) 

where {W(t)}t>o is a Brownian motion, with covariance Q, taking values in the 
Hilbert space L 2 (0, L) with the usual scalar product (■, •). We assume that the 
trace Tr(Q) is finite. So W is noise colored in space but white in time. 
Applying the Ito's formula, we obtain 

^eZ||u|| 2 = (u, dW) + [(«, u xx - uu x ) + ^a 2 Tr(Q)]dt. (43) 

as before (u, uu x ) = 0. Thus 

-^E|M| 2 = -2\\u x f + a 2 Tr(Q). (44) 

By the Poincare inequality ||n|| 2 < c||ua;|| 2 for some positive constant depending 
only on the length L, we have 

^E\\u\\ 2 < -^\\u\\ 2 + a 2 Tr(Q). (45) 

Then using the Gronwall inequality, we finally get 

E|M| 2 < E||u ||V§* + ^ca 2 Tr(Q)[l - e _ §*]. (46) 

Note that the first term in this estimate involves with initial data, and the 
second term involves with the noise intensity a as well as the trace of the noise 
covariance. 

We now consider multiplicative body noise forcing. 

dtu + u ■ d x u = ud 2 u + auw, (47) 

with the same boundary condition and initial condition as above, where wt is 
a scalar Brownian motion. So w is noise homogeneous in space but white in 
time. 

By the Ito's formula, we obtain 

1 1 

-d||n|| 2 = (u, audw) + [(u, vu xx — uu x ) + -^a 2 \\u\\ 2 ]dt. (48) 



Thus 



^E||n|| 2 = -2u\\u x \\ 2 + a 2 |M| 2 < (a 2 - ^)||u|| 2 . (49) 



Therefore, 

E\\u\\ 2 < ElluollV^ - ^*. (50) 
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Note here that the multiplicative noise affects the mean energy growth or 
decay rate, while the additive noise affects the mean energy upper bound. 
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